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We investigate the stability of the extra dimensions in a warped, codimension two braneworld that 
is based upon an Einstein-Maxwell-dilaton theory with a non-vanishing scalar field potential. The 
braneworld solution has two 3-branes, which are located at the positions of the conical singularities. 
For this type of brane solution the relative positions of the branes (the shape modulus) are deter- 
mined via the tension-deficit relations, if the brane tensions are fixed. However, the volume of the 
extra dimensions (the volume modulus) is not fixed in the context of the classical theory, implying 
\^ . we should take quantum corrections into account. Hence, we discuss the one-loop effective potential 

' of the volume modulus for a massless, minimally coupled scalar field. Given the scale invariance of 

\ the background solution, the form of the modulus effective potential can only be determined from 

■ the sign of the logarithmic term in the effective potential that depends on the renormalization scale. 
This term can be evaluated via heat kernel analysis and we show that in most cases the volume 
modulus is stabilized. In the actual evaluation, due to a lack of knowledge of the UV contributions 
from the conical branes in a six-dimensional spacetime, we consider its four-dimensional counterpart. 

' We then go on to discuss the mass scale of the modulus itself and find that it becomes comparable 
to the gravitational scale for mild degrees of warping, when the renormalization scale is set nom- 

. inally to the gravitational scale. Then, we make some suggestions on the original six-dimensional 

^ ■ model. Finally, we close this article, after discussing some phenomenological implications relating 

OO ' to the hierarchy problem of the fundamental energy scales and the smallness of the effective vacuum 

CO , energy on the brane. We find that one-loop corrections in this model appear to alleviate some of 

■ these problems. 
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I. INTRODUCTION 



CD 

String theory suggests that our universe is not actuaUy four-dimensional, but in fact a submanifold (brane) embedded 
.!Zh into a higher-dimensional spacetime (bulk). Braneworld gravity and cosmology, especially based on the proposal by 
Randall and Sundrum has developed a lot in the literature in the codimension one context String theory; 

?H however, suggests that there are as many as six or seven extra dimensions and thus, one may consider braneworld 
. . . i models with higher codimensions. Recently, codimension two braneworld has been investigated eagerly, because it 
may give a more natural resolution of the cosmological constant problem than the codimension one case does, namely 
the vacuum energy of the brane may affect only the extra dimensions, not the geometry on the brane, see e.g., 
0, 0, 0, 0, S 0> llCt frUjlj. Il3| . The basic nature of gravity and cosmology in codimension two braneworld has been 
reviewed, e.g., in Ref. [lj| (see also references therein). 

The stability of the extra dimensions is one of the most significant issues in braneworld models with two branes. 
The interbrane distance between the branes appears as a scalar degree of freedom in the four-dimensional effective 
theory and affects the geometry and cosmology on the brane. For instance, in RS-type codimension one brane models, 
this modulus, called the radion, behaves as a scalar degree of freedom in an effective Brans-Dicke gravity [l5l |. Several 
stabilization mechanisms have been discussed in the RS model, both by classical dynamics [la . IitI and quantum 
corrections to the bulk vacuum state [l9l |20| . Stability of the extra dimensions is also a significant issue in higher 
codimensions. In this paper, we focus on the issue of quantum mechanical stabilization of the modulus in codimension 
two braneworld. 
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Various codimension two braneworld solutions with two branes have been found, where the compact extra dimen- 
sions are supported by the dynamics of flux fields (4 idi B 0. • Solutions with warped compact extra dimensions 
have also been found 0, |23, 113, 13 , especially in the context of six-dimensional supergravity theory j2^ I2E ■ 
In this paper, we focus on the solution discussed by Aghababaie et al. in ^251. Similar warped solutions have been 
recently discussed in [28ll29j. in the context of six-dimensional Einstein-Maxwell theory. In these two-brane models, 
the magnetic flux plays an essential role in order to obtain regular warped solutions. By "regular" we mean that 
singularities which are stronger than conical ones are not permitted. The branes are situated at the conical defects, 
whose tensions are determined by the conical deficit angles. In warped brane models the conical defects correspond 
to the horizons of the bulk geometry and thus, the deficit angles depend on the parameters of the bulk geometry, 
e.g., the mass and charge. By the tension-deficit relations, if the brane tensions are fixed, then part of bulk geometry 
is also fixed. In the models discussed in [l^ll^, the bulk geometry is also completely determined. However, in the 
model we shall consider here |22| , which is based on six-dimensional supergravity, only the warp factor is fixed by the 
classical analysis. Thus, to fix the size of the extra dimensions, quantum corrections of the bulk field should be taken 
into account. In this paper, we consider such a supersymmetric warped fiux compactification model. ^ 

We shall consider the perturbations of a massless scalar field on this brane background and calculate the effective 
potential of the volume modulus, focusing on whether or not one-loop quantum corrections can stabilize the absolute 
size of the extra dimensions. Due to the scale invariance of the background field theory, the form of the modulus 
effective potential is almost completely fixed and therefore we only need to evaluate the sign of the coefficient of 
the logarithmic term in the effective potential to investigate the possibility of volume stabilization. This is directly 
related to the appropriate heat kernel coefficient, which is composed of contributions from both the bulk and the 
conical branes. 

The paper is organized as follows: In Sec. II, we present a warped fiux compactification model that is based on an 
Einstein-Maxwell-dilaton theory with a non-vanishing scalar field potential. In Sec. Ill, we consider the perturbation 
of a massless, minimally coupled bulk scalar field in order to derive the modulus effective potential and then investigate 
whether or not the effective potential has a minimum. As a result, the stability of the modulus is determined by the 
appropriate heat kernel coefficient. Then, we present results for our analysis of the heat kernel coefficients and we 
also discuss the rigidity of the stabilization, i.e., is the modulus mass stable against KK perturbations. In the actual 
investigation of stability, rather than the original six-dimensional model, we consider a four-dimcnsional counterpart, 
because of the difficulties associated with evaluating the UV contributions from the conical branes. Then, we make 
some suggestions for the original six-dimensional model. In Sec V, we discuss phenomenological implications relating 
to the hierarchy of the fundamental energy scales and the effective vacuum energy realized on the brane. In Sec. VI, 
we summarize this paper and discuss possible extensions of our present analyses. In Appendix A, we introduce the 
four-dimcnsional version of the warped flux compactification solution and consider the nature of the massless scalar 
field perturbations on such a background, including quantum effects. In Appendix B, we show the conformal invariance 
of the a4(/ = 1) heat kernel coefficient. In Appendix C we use the conformal invariance of the a^^coneif = 1) heat 
kernel coefficient on the cone to find the smeared coefficient a4^cone(/), which is required for the cocycle function. In 
Appendix D we apply the WKB method to give an estimate for the zeta function in four dimensions, which contributes 
to the mass scale of the modulus. 



II. A WARPED CODIMENSION TWO BRANE MODEL WITH FLUX COMPACTIFICATION 



A. Solution 



We consider a six-dimensional Einstein-Maxwell-dilaton theory with a non- vanishing scalar potential |33| as 

Se - Mi I d'x^ Qi? - ia.v^aV - le-^F,jF^' - 25?e'^) , (2.1) 

where is a dilaton field, Fij represents a U{1) gauge field strength and gi represents the dilaton potential. This 
theory corresponds to the bosonic part of the Salam-Sezgin, six-dimensional gauged supergravity theory '25', '26] with 
vanishing Kalb-Ramond 2-form and chiral scalar fields. Hereafter we set M| = 1 for simplicity and if needed, we put 
it back in explicitly. 



Recent ly, scalar and tensor perturbations about supersymmetric codimension two brane models have been discussed, see e.g., l24l 131 
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The theory given in Eq. (|2.1|l has a dilatonic, charged black hole solution 'ss'l with metric 

ds^ = -h{p)dt^ + -rj— + {2p){dx\ + dxl + dxl + dxl) , 
fi[p) 

^, , 2M i2 2 
ifiip) = -H2p), 

F,, = ^ , (2.2) 

where M and ^4 are integration constants. 

At this stage it will be convenient to perform a double Wick rotation, t — > i9, xi — > —ir with the reparameterizations 
M -M, A-^iA^. This leads to 

do'^ 

ds"^ = h(p)d0'^ + + (2p)(~dT'^ + dxl + dxl + dxl) , 

Hp) 

2M A\ ^ g\p 
p p-' 
V{p) = - ln(2p) , 

Fep^-j,. (2.3) 
In this scenario, the branes are located at positions which are determined by the horizon condition h{p) = 0; 



gl \ V 4M2 



Pi = ^(l±\/1-^1 ■ (2-4) 



We can then rewrite 



Kp)-^AP\-P')ip'-P'-) (2-5) 
and whence we also obtain the following useful relations 

+ - — , p\pl = — . 2.6 
91 91 

These give a direct relation between the mass and charge of the magnetic flux. The global period of 6 is called A6', 

which is determined by the brane tensions through the tension-deficit relations as will be discussed later. Note that 
for the p_ — > limit, we obtain a single brane solution with a naked curvature singularity in the bulk. 
By changing coordinates to 



/ 2 2\ 1/2 



Eq. (|2.3|) becomes 

ds"^ = 2p 



, 5?(pl-p^-)%2,,2^^^^,^,,^. 



(2.7) 



_.gf(l + z2)2 4(p^ + piz2) 
Thus, there are two cones at z ^ 0, C50, whose deficit angles are given by 

2 2 1 2 

5+ := 27r- ^(l-r2)A6', 5_ 2tt - ^^^M , (2.9) 

2 2 

respectively. Once the brane tensions are fixed, then r ^ p_/ p^ is fixed. But the absolute size of the extra dimensions, 
i.e., p+, is not fixed due to scale invariance of the bulk solution, unless quantum corrections are taken into account. 
In the following, instead of fixing a±, we shall fix 5+ and r. 
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B. Bulk geometry and branes 



The branes are codimension two boundaries, embedded onto conical deficits at p — p±, whose actions are given by 

S± = - J d'^xVha± , (2.10) 

respectively, where a± denotes the brane tensions. Note that we do not assume any coupling of the dilaton with the 
brane tensions and such vanishing dilaton couplings are required in order that effective cosmological constant on the 
branes vanishes completely, at least at the classical level |22] . We also assume that there is no ordinary matter on the 
branes other than their tensions. The brane tensions are related to the conical deficit angles by 

a± - MiS± . (2.11) 

We stress that these relations are only valid for sufficiently small brane tensions in comparison with the bulk scale 
Mq. Then, the angle period Ad is determined by 

, „ 27r - (5+ 27r - 5_ , , 

A9 = - = , (2.12) 

where k± denotes the surface gravities on each brane: 



^±=Tlh'{p±)^^{pl-p^_). (2.13) 



Thus, we obtain 



After eliminating r, we obtain 



27r-(54 
27r-(5_ 



(2.14) 



AO.'-^^^^li^. (2.15) 

We emphasize that, in general, A6 ^ Itx implying that the bulk angular period can be completely determined by the 
brane tensions. For later convenience, we regard it as a function of r and 5+: 

A%,^+) = ^^P^. (2.16) 
3^(1 -H) 

Once the brane tensions, tT+ and U- are fixed, then r is also fixed and thus, we now regard the free parameters 
as r and 5+ , along with the dilaton bulk coupling g\ . The remaining degree of freedom used to determine the bulk 
geometry is the absolute size of the bulk, i.e., p+. Due to the scale invariance of the background brane solution, 
can only be fixed by quantum corrections of bulk fields. Note that there is also magnetic flux constraint given by 

dp r deF,e = -AM(^ - ^) = , (2.17) 

P- Jo P- P+ gir 

but the magnetic flux only depends on r and Sj^ and docs not contribute to the volume of the bulk. 



C. The volume modulus 



We consider the modulus dynamics in a four-dimensional model in the context of the moduli approximation, as- 
suming p+ is a function of x^^. This approximation may be valid for sufficiently low energy scales compared with 
the bulk gravitational one Afg. Like for the codimension one case, we expect that this approximation is valid for low 
energies, namely when the energy scale on the brane is below that of the bulk. In this subsection, we shall reintroduce 
Mq in the action in order to clarify the dimensionality. We regard p-f. as a function of the coordinate on the brane 
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= (r, X2,X3, X4) in Eq. H2.2|l and thus, as the volume modulus. Then, we expand the components of the Lagrangian 
density and take the modulus parts out. It is rather useful to move to the four-dimensional physical frame: 



g2(l + z2)2 4(1+^.2^2)2 

2p+(x^)p(z) 



dz^ , 5?(l-r2)2 



= 2p{z) p+{i^^) 



52(l + z2)2 4(1 + ^.2^2)2 



+ 



.g2(l + z2)2 4(l + r2z2) 



r/fj^i^dx^dx'^ 
z-d0'^ + -q^udx^dx" 

+ rj^^dx'^'^dx^ 



z^de^ 



(2.18) 



where has the dimension of physical length. The metric form is just like a Kaluza-Klein theory with two compact 
dimensions, apart from the overall conformal (warp) factor 2p{z). 
Then, we evaluate the effective action of the volume modulus 



-g{R-{d^f) 



(l-r2)z 

mod - p^(l + ^2)2 



{dp. 



4p+d'p+ 



(2.19) 



Note that the flux and dilaton potential terms do not give any kinetic term contributions to the moduli kinetic term 
because the vector potential Aj has no scale dependence. The second term is a total derivative term and does not 
contribute to the effective action. Integrating over the extra dimensions, we obtain 



Af|(27r 



29i 



d^x 



[dp+f 
P+ 



After redefining the modulus as 



X6(2:n 



/4M4(27r-5+) 



9l 



P+ 



we obtain the canonical form of the modulus kinetic term as 



mod 



= / d^x 



(4(^Xa) 



(2.20) 



(2.21) 



(2.22) 



III. THE ONE-LOOP EFFECTIVE POTENTIAL OF THE VOLUME MODULUS 



Next, we introduce a massless, minimally coupled scalar field (not to be confused with the dilaton field that supports 
the background configuration). From this we can investigate the one-loop effective action for such a scalar field and 
thus, the effective potential of the volume modulus on the co-dimension two warped brane background. From now 
on, we shall work mainly in Euclidean space, rather than in the original Lorentzian frame. 

The action for the massless scalar field perturbations is given by 



^scalar 



4 / '^^^Vff'^^6<; 



(3.1) 



A. Scalar one-loop effective action 



The one-loop effective action for a massless minimally coupled scalar field is defined as 

We = ilndct(-A6), 



(3.2) 



where Ag is the six-dimensional Laplacian, which is divergent and needs to be regularized and renormalized. For this 
purpose, we define 



u2« /■°° dt 



(3.3) 
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where K{t, Ag) is the (integrated) heat kernel defined by 

i^(<,A6) = Tr(e-*(-'^''')) . (3.4) 
The (integrated) zeta function is related to the heat kernel by a Mellin transform: 

C(s, As) - dtt^-'K{t, Ae) = Tr((-A6)^) (3.5) 

and after analytically continuing to s ^ we obtain the renormalized one-loop effective action. The explicit expression 
for the renormalized effective action is 

- -^^i'^T{s)as, As) = -i Q - 7 + In/i^^ C(0, Ag) - ^('(0, Ag) + 0{s) , (3.6) 

where 7 is Euler's constant 7 « 0.577216 and the pole term at s = is removed by renormalization. Thus, after a 
redefinition of the renormalization scale we obtain the renormalized effective action and the renormalized scalar field 
effective action can be written as 

We,rcn = -ic'(0, Ae) - ^((0, Ag) In^^ _ (3 7) 

By integrating over the internal dimensions, the four-dimensional effective potential is 

W^6,rcn= J {d''xpl)V6,eS = J d^iV^^oS, (3.8) 

where Vcs has the dimensions (length)^'^. In the next section, we shall derive the effective potential in a codimcnsion 
two warped brane model. For brevity, from now on we shall omit the subscript "ren" . 
In the case of codimension two, the zeta function is given by the summation 



d*k 1 



)4 A2 



(3.9) 



where the eigenvalues are defined by 
It is straightforward to show that 



AecPx^-y^bx. (3.10) 



C(0,A6)=a6(/ = l), (3.11) 
where aQ{f ) is a heat kernel coefficient, defined by the asymptotic expansion of the heat kernel [s^. l35l|: 

if(i,A6)c.^i('=-«)/2a,(/), t^O. (3.12) 

k>0 

In the warped codimension two case, the eigenvalues are unknown and there are not even any analytic solutions for 
the eigenfunctions so we must resort to using an approximate WKB method to estimate the effective action. 

B. Continuous conformal transformations 

One strategy to evaluate the one-loop effective action and the effective potential is to define a continuous conformal 
transformation (parameterized by e) 

dsl^ = e^i^-'^-dsl, c. = iln(2p) (3.13) 

and thus 

d~4 = i^PT ( 2^2 + ^IK^'f^lf^ dO' + d^A , (3.14) 
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where for e = 1 we have the original metric, which we shall denote as Ag^e = Ag. The classical action of this scalar 
field is changed under a conformal transformation, see Eq. (|3.13|) . by (here we are considering a massless, minimally 
coupled bulk scalar field) 

5scalar=-^y" <f X^4>^e4> = j (Ag + (e)) , (3.15) 



where 



EQ{e) = -4(e- l)2g'^''Vat^VhCj + 2(e- l)A6lncj 

1 \ ^ g\l - e){p\ p^_){pI{2 + (1 - e)z^) + p'^.z^-l + e 2z^)] 



2p) {pl+plz^y 



(3.16) 



Due to technical reasons, which will be explained later, we shall evaluate the zeta function in the non-warped frame 
e = 0. The correction associated with such a conformal transformation is commonly known as the cocycle function 
(obtained from integration along the paths of the conformal transformation): 

We - -ic'(0,A6)-ic(0,A6)lnAi2 = -ic'(0,A6,,=o)-^C(0,A6^.=o)ln/i'-^ deae (/ = 9e In^) ■ (3.17) 
Furthermore, thanks to the conformal invariance of C(0, Ag) = aQ{f — 1) we arrive at 

We = -ic'(0,A6,=o)-^C(0,A6,,=o)ln/i'-^ de {f ^ d,{e - l)u;) 

= -ic'(0,A6,,=o)-ia6(/ = l)lnM'-^ de ag (/ = i ln(2p)) 

= -laeif = 1) Hp'p+) + { - ic'(0, A6,.=o) - de {ae{f = \ ln(2p)) - 1 ln(p+)a6(/ = 1)) } 

= -ia6(/ = l)ln(MV) + {-^C'(0,A6,,=o)-^'*a6(/ = ^ln(|^))}. (3.18) 



The term ag is given by the volume integration of linear combinations of cubic order curvature tensors | 34l l35l l36| : 

p Tyj pkn I p p Tjikjl p pj pknli , p pij ryknlp , p pi fe TDjlnp\ 

-—R^.R^R +-R,RuR ^-R,kRnuR +^R,knRipR + ^R.knR i pR } 

oDU 



-mElR + UEeR'^.k + ^EeR^ - 2EeR^ R^. + 2EeR'^^R^ 



i]kij 



-(contribution of conical branes)| . (3.19) 



C. Effective potential of the volume modulus 

In order to avoid a volume divergence we shall employ the one-loop effective potential, as given by Eq. (|3.8|) . to 
try and stabilize the extra dimensions. From our discussion in the previous subsection, the effective potential takes 
exclusively the following form 

T. , r X A6(r, S+) - Bejr, S+) \n{p^p+) 

V6,cff(r,d+,p+,M) = 2 ' (3-20) 

P+ 
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where we have defined 

d'xAeir, 5+) = f d^jM^ ^ _ r dea^ - i \n{^)) - ^CiO, A6,.=o) , 
J P+ Jo 2 p+ 2 

J d'xBeir, S+) = J d^jM^ ^ i^(o, Ae^e^o) = ^aeif = 1) . (3.21) 
Clearly, if Bg^r^ (5+) > 0, then the modulus effective potential has a minimum at 

PI = ^-2e(2Aa + B6)/(2B6) (3 22) 

After a redefinition of the modulus, as given by Eq. (|2.21|) . the effective potential can be rewritten as 



' 4Mi{27r - 6+) \^Mr, S+) - B^jr, 5+) In ( ^"^'^g 
^ 9^ 

The field value at the minimum is then given by 



V6,cff(f,d+,x,M) ( 2 4 ■ (^-^^^ 



^ 4M|(2^ -,5+) ^(2A.+B.)/(2i^.) (3 24) 



and hence, the effective modulus mass becomes 



■■= 5L^6,cff[ = ^^i?6(r,<5+)e-^(^^«+^«'/(^^«)(^)'A/| . (3.25) 



Thus, the ratio between the stabilized mass and the effective bulk cosmological constant is 



.B6(r,5+)f-^^)V(2^«+^^'/^«. (3.26) 



mjcs ^ 8_ 

The expression above allows us to quantify the rigidity of stabilization; namely, is the modulus mass stable against 
KK perturbations. However, the problem is that the og contribution for conical branes has never been formulated as 
far as the authors are aware. 

IV. VOLUME STABILIZATION 

In this section, in order investigate the one-loop effective potential of the volume modulus, we consider the four- 
dimensional counterpart of the six-dimensional case in the theory |33j| . 

^4 = Mi J d^x^g {r - ia,^9V - ^e-^i^,,i^*J" - Ag^'P^ , (4.1) 

where, again, is a dilaton, Fij represents a U{1) gauge field strength and g represents the dilaton potential (cor- 
responding to gi in the original model). Hereafter, we set the bulk gravitational scale to M| = 1 for simplicity; we 
shall reinsert it if and when it is needed. The essential properties of this solution is summarized in Appendix A. The 
reason we use this model is that we can evaluate the UV contributions from the conical branes. We shall extrapolate 
our results to the original model in six dimensions. 

A. One-loop effective potential of the volume modulus 

The zeta function and one-loop effective action in four dimensions can be defined similarly to the case of six 
dimensions, as in Eq. H3.9I) and Eq. (|3.8|l . just by replacing "6" — > "4" (and (27r)^ —> (27r)^). Thus, the one-loop 
effective potential of the volume modulus takes the same form as Eq. 13.8|l . after the replacement of d'^x with d'^x. The 
effective action also takes the same form as Eq. (|3.18|l . given the scale invariant nature of the flux compactification 
solution, implying the effective potential can be written similarly as 

,r ( , N Ai{r,5+) - Bi{r,5+)\Yi{p'^p+) 

V4,off(r,()+,p+,M) = , (4.2) 

P+ 
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where 

P+ Jo 2 p+ 2 

d^xB^{r. 5+) = f d^^MlM = 1^(0, A4.e=o) = ^a^f = 1) • (4.3) 
J P+ 2 2 

Note that the stabiUzed volume of the internal space is given by 

= ^-ie(A4+s,)/(2B,) _ (4 4) 

Again, if Bi{r, (5+) > 0, then the modulus effective potential has a minimum at p-^- = p+,*, which is determined by 
v44(r, (5+) and Bi{r,5^). Importantly, given the similarity of the spacetime structure in this four dimensional model 
with the original six dimensional model, the form of the effective modulus field can be estimated by dimensional 
arguments, in comparison to Eq. H3.24|l . Namely, 



X4(x^):^W ^'^2 —P+ (4-5) 

acts as the modulus field in two dimensions. The field value at the minimum is then given by 

and hence, the effective modulus mass is 

mU ■■= = ^^S4(r,5+)e-2(^^+^^)/^^(^)\/|. (4.7) 

The rigidity of the stability is determined by the quantity 

m 



'off 



52(p+,,)-i 211 -5, 



S4(r,^+)e-(^^+^*)/^^(^) . (4.8) 



Given the conformal invariance of the heat kernel coefficients, 04(7) is straightforward to evaluate in the e — 1 
frame. In the original four-dimensional spacetime, including the contribution of the cones |37l ISSl Is^ we have 

aiU = 1) := (47r)-2|360-i ^ d^x^i^UR'"" ,k + 5i?' - 2R,R' + 2R,,kiR'''') 

^ /^^-E (|)^^[fG^+A,5:(.H^-24"^^,^)) 

A— di 27r a 

+ ' n "^-tf ( E - 2 E ^-^-^ + \ E -'^'^ + E {-'^'^ " I (4-9) 

^ 27r / a a,b a a 

where {a} runs from 1 to 2 which denotes the two orthonormal directions to the conical branes 



and is the extrinsic curvature defined by = — ft-^/i' Vfcn|°\ Here denotes two mutually orthogonal (inward 

pointing) unit normals and a ^ z, 9 and k*-"-* — h^^ K^j^\ where hij = gij — X^a induced metric on the 

brane. The terms Ai,2 are undetermined parameters associated with conformal invariance of the conical heat kernels. 
But, as shown in Appendix C, the extrinsic curvatures vanish on the cones and we need not specify them. Note that 
this formula is valid for any deficit angle (5^ l39| . 
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FIG. 1: The bulk contribution to -84(7-, 0.01) is shown as 
a function of r, for g = 1. 



FIG. 2: Contribution of conical branes to -B4(r, 0.01) is 
shown as a function of r, for g = 1. 



B. Stability of the volume modulus in four dimensions 



Now we present the results of our analyses on the one-loop cfFective potential of the modulus. The function i?4(r, (5+) 
that corresponds to BQ{r, (5+) in the original six-dimensional model is composed of two parts; a bulk and conical brane 
contribution and thus, the total is 



(4.11) 



From the properties of the heat kernel coefficients i?4(r, (5+) is invariant under conformal transformations, as is 
explicitly demonstrated in Appendix C. The bulk part of i?4(r, S^) is obtained from Eq. ljB2|l for e — 1 and as shown 
in Appendix B, it is in fact independent of e. Similarly, the brane contribution to B4{r, 6+) is obtained from Eq. (|B8|I . 

In Fig. 1, the bulk contribution to B4{r, S^) is shown as a function of x, for g = 1 and S+ = 0.01. Note that in the 
following plots, which are described as functions of r, we have set 5+ = 0.01. There is no significant dependence on 
(5+ and thus we mainly focus on r. For other values of g, the behavior of B^lr, 6+) is also almost the same. The bulk 
contribution is always positive and as r ^ it suffers from a divergence because of the presence of a naked curvature 
singularity. In Fig. 2, the contribution from the branes, to i?4(r, ^+), is shown as a function of r. Comparing this 
with Fig. 1, one can readily see that the cone contribution becomes comparable to that of the bulk one but the sign 
is opposite. Thus, the positivity of the total, i?4(r, 5+), is a rather non-trivial issue. In Fig. 3, the total, B4{r,6+), 
is shown as a function of r and we find that the sign of i?4(r, 5+) is always positive. Thus, for all values of r, the 
modulus can be stabilized; this is one of the main results in this article. Though here we discuss the four-dimensional 
counterpart, we will make suggestions on the more realistic case of six dimensions later. Note that the behavior as 
shown in Fig. 3 is almost the same for other values of the conical deficit angle 5+ . 

Next, we wish to discuss the minimum potential energy and and thus, the effective mass of the modulus. These 
can only be considered through the combination ^4(7-, ^+)/i?4(r, 5+), sec Eq. (14. 4f) . where A4{r,6+) corresponds to 
AQ{r,d+) in the original model. This ratio should then be compared with the rcnormalization scale, /i. We give an 
explicit derivation of the cocycle function and the derivative of the zeta function in Appendix C and D, respectively. 

The term A4(r, (5+) is obtained from Eq. H4.3|) . The first term is the derivative of the zeta function in the non- warped 
conformal frame, e = 0, in the context of WKB approximation. The second term is known as the cocycle function, 
which corresponds to corrections associated with a conformal transformation. By combining these two contributions 
we obtain 



A4(r,5+) = (A4(r,5+))_^^,^+(A4(r,^+)) 



'WKB- 



(4.12) 



The bulk and brane contribution to the cocycle function, A4.cocycie('', ^-1-), is given in Eq. (jClip . while the WKB 
contribution, A4^^kb{i", in the non-warped frame, is given in Eq. ljD38|) . Actually, the value of A4{r,S+) itself is 
not as important as the ratio A4(r, 5+)/ B4{r, J+), because the stabilized size of the extra dimensions is given by Eq. 
(|4.4|) . In Fig. 4, we show a logarithmic plot of the exponential factor as a function of r for g — 1 and 5+ = 0.01. For 
all values of r, the stabilized size of the extra dimensions is larger than that of the rcnormalization scale. Especially 
in the limits r and r ^ 1, much larger extra dimensions are permitted. 

Furthermore, in order to show how the stability of the extra dimensions are set rigid, we can consider the ratio 
between the effective modulus mass and the bulk curvature scales g'^{p^^^)~^ . In Fig. 5, we show logarithmic plots 
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FIG. 3; The total logio(_B4(r, 0.01)) is shown. The red 
(sohd), green (dashed, with wider intervals) and blue 
(dashed, with shorter intervals) curves correspond to 
g = 0.5, 5, 50, respectively. 



FIG. 4: A logarithmic plot of the exponential factor in 
Eq. (14.411 is shown as a function of r for g — 1 and 
5+ = 0.01. In the following all logarithmic plots are to 
the base 10. 




FIG. 5; A logarithmic plot of Eq ijH^ for ^ = M4 is 
shown as a function of r. The red (solid), green (dashed, 
with wider intervals) and blue (dashed, with shorter in- 
tervals) curves correspond to g = 0.5, 5, 50, respectively 
and S+ = 0.01. 



FIG. 6: A logarithmic plot of Eq igSl for fi = M4 is 
shown as a function of S+ . The red (solid) , green (dashed, 
with wider intervals) and blue (dashed, with shorter in- 
tervals) curves correspond to g = 0.5, 5, 50, respectively 
and r — 0.5. 



of Eq. (|4.8|) as a function of r for various dilaton couplings g. In Fig. 6, we also show them as a function of If 
the ratio is greater than unity, then, the stability of the modulus is generic. From Fig. 5 and 6, the value of Eq. H4.8() 
becomes larger for larger couplings g (in fact the relative mass is almost proportional to g^). Thus, for large dilaton 
couplings the stability of the modulus is rigid against KK perturbations. In Fig. 5, we also see that for smaller values 
of r, the modulus mass becomes relatively light, whereas for larger r it is insensitive to r. In the limit r the 
modulus mass diverges, because a bulk naked singularity is formed. As a final comment. Fig. 6 is not particularly 
sensitive to changes in (5+, but it does exhibit a divergent behavior for (5+ 27r. This is because in such a limit the 
bulk essentially disappears. 



C. Suggestions for volume stabilization in six dimensions 

Now, we shall make some suggestions on the stability of the volume modulus in six dimensions. Because of our poor 
knowledge about the contribution from the conical branes we can evaluate only the bulk contribution. However, we 
can compare the result of the bulk contribution in six-dimensions to the four-dimensional one discussed previously. 

Evaluating the bulk part of -65(^7 5+) from Eqs. (|3.19|) and H3.21|l leads to 

(r.(.^^ ^ 1 r"" 1 gti2TT-S+)z 1 ^{l,r,z) 

l^6ir-,o+jjj^^j^ 2 io 64^2 2(1 + ^2)(i+ ^2^2) 1260 (1 + r2z2)6 ' ^^'^^^^ 
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FIG. 7: The bulk part of Be (r, 0.01) is shown as a func- 
tion of r for gi = 1. 



FIG. 8: logjQ(7?(r, 0.01)) is shown as a function of r, 
for 5+ — 0.01. The red (solid), green (dashed, with wider 
intervals) and blue (dashed, with shorter intervals) curves 
correspond to g — 0.5, 5, 50, respectively. 



where 



*(P+,P-,^) 



2035 p^^z^ + 1564 pip\_ - 548 p^_p\ 
5461 p^^z^ - 1740 p^"p- + 4964 p 



5461 pL^z* + 4964 pi^z'*' + 2035 pi^z'' + 1748 pi^z 



^12.^10 



12^6 



^12.^12 



12^2 



12757 p^z>* - 2882 p^'z>l - 9049 p+'^zVl - 5936 p+p-^^ - 548 p^z'>^ 
1740 z^VIp-" + 711 pV^^V- + 3652 p^z'^^pi - 9049 z V+P-° - 2882 p"'p+^' 



1748 p^^ + 32534 p>_z* + 26109 p+zV* + 13944 p^p 

1564: z^^plpt 
nip^^pW 



5936 zV^pl" 



32534 z^p^p^ 

„6 „6 ,6 



13944 zi°ptP- 



26109 p^p^z^ 



3652 plp+z^ + 12757 p^.p^z^ 



,6 „6 



10480 p^p^z^ - 30044 p^p^z" - 27054 p^z^^ - 10480 p>^z"' - 27054 p^p' 



(4.14) 



In Fig. 7, we plot Eq. (|4.13|) as a function of r for gi — \ and — 0.01 (this behavior is essentially the same for 
other values of the coupling 51). The behavior of the bulk part of BQ{r, 5+) is similar to that of Bi{r^ 5+) as shown 
in Fig. 1. As we mentioned previously, we have no analytic formulae to evaluate the contribution from the conical 
branes for ae{f). However, by analogy we expect similar behavior for the conical part of the heat kernel coefficient. 
It may be worth mentioning that for larger r and appropriate values of gi, BQ[r, (5+) is very likely to be positive and 
thus the modulus can probably be stabilized. We could also evaluate C'(0, Ag), via the WKB method, and due to 
symmetry would follow almost identical steps to the four dimensional case, see Appendix C. However, we have no 
way to include the cocycle correction, at the level of our semi-analytic approach. 

Thus, any further investigation of the conical contribution and the effective mass of the stabilized modulus (in six 
dimensions) is somewhat out of the scope of this article and is left for future work. 

V. PHENOMENOLOGICAL IMPLICATIONS 

Finally, before closing this article, we shall discuss some phcnomcnological implications, i.e., the hierarchy struc- 
ture between the fundamental energy scales and the effective vacuum energy on the brane. Even though we are 
investigating, primarily, the four-dimensional case, we expect similar results to hold in six dimensions. 



A. Hierarchy of the fundamental scales 



One of the most longstanding issues in theoretical physics to date, is the problem of the hierarchy structure of the 
fundamental energy scales, i.e., the huge gap between the gravitational scale, Mpi ~ 10-^^GeV, and the electroweak 
scale. Mew ^ lO^GeV. 
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In the original six-dimensional model, the effective four-dimensional Planck scale^ is 

Ml,^P±^^^^Ml. (5.1) 

If we assume a brane localized field whose bare mass is given by m? on either brane at p± then the observed mass 
scales are 

2 2 2 2 2 Ir 

mj_ = m , m_ — r m . (5.2) 

Thus, the mass ratio between the field and the effective Planck mass is given by 

m|_ _ ( l^]f^\_^_ -(2Ae + Be)/(2Be) !!f^ ^ f [^Jff_\2 ^l_-{2Ae+Be)/(2Be) (r: o\ 

~ ^ ^6 /27r-,5+ ' - V A/4 r 2TT-S+ ' ^ ^ 

Assuming that the factor of {jim/MQ)'^ takes the optimal value of 0(1) at the unification of the fundamental scales, 
the mass ratio becomes 

_3l_g-(2A6+B6)/(2B6) ^ (5 4) 

27r — (5+ ' 

where we have used the value of given by Eq. (|3.22(l . at stabilization. At present, the best we can do is use the 
results obtained from the four-dimensional model: 

R(^r, 5+) = ^^e'^A,+B,)/B, (5 5) 

27r — 0+ 

The results are shown in Fig. 8 and indicate that smaller values of g and r lead to a large hierarchy on both branes. 
However, the mass of the stabilized modulus is perhaps too light to support the realization of the hierarchy at the 
fundamental scales. This is somewhat similar to the case of RS codimension one branes Il9ll20ll. 



B. Effective vacuum energy density 

Another interesting issue is the origin of the observed vacuum energy density, which is much smaller than the field 
theoretical prediction; the so-called the cosmological constant problem. Given the background brane solution that we 
are considering, the brane vacuum energy (the brane tensions) only contributes to the bulk geometry and therefore 
does not act as an energy source for today's accelerated cosmic expansion (at the background level the branes are to be 
kept flat). However, in this model, the observed vacuum energy density can arise from one- loop quantum corrections 
as follows. 

At the stabilized volume, from Eqs. (|3.2U|) and l|3.22|l . the effective potential takes the value 

Vl^^ = -i/S6(r, 5^)e-(2^a+^a)/(2^^) (5.6) 

and hence, the brane realized vacuum energy is almost completely determined by the renormalization scale. It is then 
quite natural that the renormalization scale is taken to the unification scale, i.e., of order Mg. However, in such a case 
the exponential factor is also very small and therefore we may be able to obtain a sufficiently small vacuum energy, 
which acts as a dark energy source. Although the vacuum energy is nega tive definite one might expect some kind of 
uplifting mechanism, for instance, due to the presence of test branes |40l |. 

As in the case of the hierarchy problem, we shall consider the potential energy at the time of the stabilization. 
Whence, the vacuum energy is given by 

- ~^Ji^B^{r, 5+)e~(^^+^-)M^ . (5.7) 

The coefficient acting on the renormalization scale is very similar to the quantity plotted in Fig. 5. Again, for smaller 
g and r, we can expect smaller vacuum energy densities of the volume modulus. The mass of the modulus is also very 
light, but this might be a sign of a dynamically varying vacuum energy. 



^ Not to be confused with M4 in tlie four-dimensional model. 
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VI. SUMMARY AND DISCUSSION 

We have discussed the stabihty of the volume of a warped codimension two brane model with flux compactification, 
by taking the one-loop quantum corrections of a massless scalar field into account. We have also examined the rigidity 
of stability by comparing the ratio between comparison the effective modulus mass and the bulk curvature scales. 

First, we introduced a warped flux compactification model in six dimensions, where we started from the solution 
with a spacetime structure similar to the Reisner-Nortstrom black hole. After a double Wick rotation, we were able 
to derive the desired braneworld solution. The branes are located at the horizon positions and correspond to conical 
singularities. For relatively low energies, namely when the brane tensions are less than the bulk gravitational energy 
scale, the brane tensions are related to the deficit angles directly. Once the brane tensions are fixed, via Eq. H2.14|l . 
the ratio of the brane positions, r = p-/ p+ is also fixed. But, the volume of the extra dimensions is not completely 
determined, because of the scale invariance of the solution. To this end, we investigated whether or not this dynamical 
degree of freedom can be fixed by considering the one-loop quantum effective action for a massless, minimally coupled 
bulk scalar field. We investigated the stability of the modulus for given values of the two brane tensions, or equivalently 
of r and 5+ (the deficit angle of the brane located at 

The form of the effective potential can be fixed up to overall coefficients, due to scale invariance and dimensional 
arguments, which depend on r and S+ as in Eq. H3.20|l . From this the condition for the volume modulus to have 
a minimum (and hence be stabilized) is determined only by the sign of the coefficient in the logarithmic term, 
i.e., Be{r,S+). This coefficient is directly related to the zeta function ^(0, Ag) and thus, the heat kernel coefficient 
a6(/ = l). 

In the actual investigation of stability, we have used an alternative model in four dimensions that has a very similar 
spacetime structure to the original six-dimensional model. This is because the contribution for the cone to the og heat 
kernel coefficient has never been formulated as far as the authors are aware. We then showed that the heat kernel 
coefficient is positive definite for most dilaton couplings, g, independently of the choice of r and 6^. The contribution 
from the bulk is positive, whereas that from the conical branes is negative for smaller r; though they are of the same 
order. For all of the cases with g ~ 0(1), the bulk dominates over the brane parts and thus, the volume is stabilized. 

However, one might ask whether or not the stability is rigid against the various KK perturbations. This rigidity 
can be determined from the ratio between the effective modulus mass and the stabilized curvature radius of the bulk, 
as given by Eq. H4.8|l . implying a knowledge of the ratio of Ae{r, S+)/BQ{r, S+). If this ratio is less than unity, then 
the mass of the volume modulus becomes too light and may be destabilized by other KK contributions, whereas if it 
is greater than unity, the modulus field is not easily perturbed from the "stabilized" minimum. We showed that the 
mass ratio becomes larger for larger coupling which is proportional to g^ and therefore, stability is rigid. For larger 
degrees of warping, i.e., for smaller value of r, other than r — 0, we found that the stabilized modulus mass becomes 
relatively light and that this ratio is not so sensitive to different values of deficit angle d+. Note that it only exhibits 
a divergence at (5+ = 2tt, where the bulk essentially vanishes. 

Then, we made some suggestions for the six-dimensional case and by comparing the bulk part with that of the 
four-dimensional 04 bulk part we found similar behavior indicating that our conclusion in four dimensions will carry 
through to the six-dimensional case: The modulus can be stabilized and furthermore our approximate WKB approach 
suggests that the mass of the modulus leads to rigidity for certain couplings and parameters. Of course, once formal 
expressions for the six-dimensional conical heat kernel coefficients are found, then we can extend the same analysis 
as that employed in this paper. However, a more accurate (semi-analytical/numerical) method to evaluate ^'(0, Ag) 
should be found, rather than the approximate WKB eigenvalue method that we used in this article. 

We also discussed some phenomenological issues, i.e., the hierarchy problem of the fundamental energy scales and 
the extremely tiny value of the observed vacuum energy density. The six-dimensional results were extrapolated 
from the corresponding results in four dimensions. We found that in the hierarchy problem, the ratio between the 
energy scale of the brane localized field and the gravitational scale becomes much smaller than unity, especially for 
larger degrees of warping r <C 1. Furthermore, the effective vacuum energy density can also be much smaller than 
the gravitational scale. The vacuum energy is negative by definition and we therefore need some kind of uplifting 
mechanism, to obtain a more realistic cosmology. In these cases, as was mentioned previously, the mass of the volume 
modulus may also be small and thus, might be destabilized. 

It is also interesting to compare our results with the codimension one case. An important difference is that in the 
codimension one case all the moduli are not dynamically fixed and should be stabilized by quantum corrections of bulk 
fields; however in the codimension two case at least one of the moduli fields is fixed dynamically, because the brane 
tensions are directly related to the bulk geometry. In terms of heat kernel coefficients another obvious difference is 
that we are considering even dimensions. In odd dimensions only boundary terms occur and thus, the renormalization 
only affects the brane tensions, whereas in even dimensions there is a bulk contribution and hence logarithmic terms 
appear. Probably the most important difference is that the stabilization of the volume modulus can be realized by 
the contribution from the bulk, not from the branes, for small warpings. For r <C 1 (large warpings) the stabilized 
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modulus effectively behaves as a very light scalar field on the brane, which is somewhat similar to the codimension 
one case, e.g., see [l9ll20j . 

Before closing this article we shall refer to some extensions of this work: One issue is to do with the quantum 
corrections of other fields, such as the magnetic fiux and graviton. The tensor part of the graviton perturbations 
should be equivalent to a minimally coupled scalar field, which we investigated in this article and such analyses 
might give us important insights into warped fiux compactifications with self-gravitating branes. Besides quantum 
efi^ects, there are many interesting issues that still remain in codimension two braneworld, e.g., gravitational waves 
in the bulk j4lj| . non- linear dynamics of the bulk and branes, thick branes and so on. Such work might be key to 
understanding the self-tuning mechanism of the effective cosmological constant as a dynamical process. We hope to 
report on these topics in ongoing work • 
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APPENDIX A: A FOUR-DIMENSIONAL COUNTERPART 



1. The model 



In this Appendix, we consider an alternative model of a warped flux compactification in a four-dimensional Einstein- 
Maxwell-dilaton theory given by Eq. (|4.1|l . This is a toy- model analogue of the original six-dimensional co-dimension 
two braneworld model with a warped fiux compactification and a non- vanishing scalar potential 33|: Hereafter, we 
set the bulk gravitational scale to Mf — 1 for simplicity; we shall reinsert it if and when it is needed. 

As for six dimensions, we first consider a black hole solution of the theory [23|- After a double Wick rotation, we 
obtain 

ds" = h{p)de'' + + {2p){~dT^ + dxl) , 
n2 

h{p) = 2M - — - 2g^p, 
P 

^(p) = -ln(2p), 

F.p^-^. (Al) 
In this scenario the branes are located at positions that are determined by the horizon condition — 0; 



2g2 \^ V M2 

We can then rewrite 



P 

and whence we also obtain the following useful relations 



Mp) = ^(p+-p)(p-P-) (A3) 



P+ + P- = ^ , P+P- = TTT • ^4 



This gives a direct relation between the mass and charge to the magnetic fiux. The global period of 0, denoted A6', 
is determined by the brane tensions through the tension-deficit relations as is discussed later. 
After a coordinate transformation 



X 1/2 

P- P- 



(A5) 
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the bulk metric becomes 

- 2p f + g' ^P^- P-^l'l de' - dr' + d.l) , (A6) 

where our choice of factorization is for later convenience. In this frame we can now easily investigate the asymptotic 
behavior of the spacetime near each of the conical branes. For z — > 0,oo, there are cones whose deficit angles are 
obtained from 

5+:=2TT-g^{l^r)Ae, S- := 2n ^ g^^—^A9 . (A7) 

r 

The modulus corresponds to the proper length between two branes (strings), which is characterized by p±, or 
equivalently p+ and r. The string tensions are related to the conical deficits by 

a± = A/|(5± . (A8) 

We stress that these relations are only valid for sufficiently small brane tensions in comparison with the bulk scale 
M4 . Again, A9 is given by a relation identical to Eq. (|2.12(l and thus. 

After eliminating r, we obtain 

Once the brane tensions, ct+ and tT_ are fixed, then r is also fixed and so we now regard the free parameters as r and 
cr+, along with the dilaton bulk coupling g. The remaining degree of freedom used to determine the bulk geometry 
is the absolute size of the bulk, i.e., p+. However, because of the scale invariance, p+ can only be fixed by quantum 
corrections of bulk fields. The magnetic flux is given by 

dp r dOFe, = -AeQ{— -—)= -^!^./?. (All) 

Thus, the gauge field flux does not contribute to the size modulus. 



2. Massless scalar field perturbations in a four-dimensional model 

We discuss the basic properties of a massless, minimally coupled scalar fleld on this background, which is introduced 
in order to investigate the one-loop quantum effects in the counterpart model discussed in Sec IV and V, whose action 
is given by 



■'scalar " 2 



i J d''x^<t>Ai^. (A12) 
As before, we consider a continuous conformal transformation of the metric, parameterized by e; 

d~sl ^ = e2(^-i)"ds^ , uo^\ ln(2p) , (A13) 



2 

and thus 



ds 

where for e = 1 we have the original metric, which we shall define as A4.e = A4. The classical action changes as 



d^x^^A^^ '^^^V^^ + Ei{e)) , (A15) 



*^scalar 

where 

Ei{e) = -(e- l)25'^''VatjVbCJ + (e- l)A4lntj 

1 V .9^(1 - - P-){p+(2 + (1 - e)z^) + P-z\-l + e - 2z2)} 



2p) {p++P-z^? 



(A16) 
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APPENDIX B: CONFORMAL INVARIANCE OF HEAT KERNEL COEFFICIENTS 

The classical action of a massless, minimally coupled scalar field is changed under this conformal transformation, 
see Eq. ()A15|) . In four dimensions, including conical branes, the 04 heat kernel coefficient is given by 

a4(/ = l) := (47r)-2j 360"M d'^Xy^ ,k + QORE4 + ISOEj + 12R''.k + - ^R^jR'^ + 2R^jkiR"^''^)} 

[ Jm 

■' A=± ^ 2n a 

- Wo Vi^itf ( E 2 ^ + I E + A. - ^4'^^ I • (Bl) 

^ 27T ' a a,h a a ) 



1. Contribution of the bulk 

The integrand of the bulk part of Eq. ljBl(l is 

V 4V , +;;buik I4407r2 7o (l + z2)(i + rz2) (l + rz2)4 ^ 

where 

*(p+,p_,e,z) = -216ezV+P- -40pi^e2zV- + 504zV^/92_ - 144 zV+ P- + 48 z**ep^p?. + 8 z*^eV+P- 

- 48 p+e^zV^ - 112 p%e^z^p^ + 192 p^p+e^z* - 96 p+e^zV?_ - 88 e^z^p^_p+ 

+ 64 p+ e'''zV^ + 33 p\z-^ + 8 + 33 piz'' + 48 ep\pi - 36 p+zV^ + 156 p+ z^- 

+ 72 ez^pt - 144 ep^p_ + 24 (^z^p\ - 16 p^e^z" + 16 e'z'^p^_p+ - 48 e^plz^pl 

+ 48 e^p5|_z2p- + 104 z'^eV+yO- " 16 z*^eV+ pt - 16 e^z^i - 16 e^z^pi - 16 p_ e^p^ 

+ 168 ezVi + 8 e^z^pt + 8 eV^p+ + 96 ez^pt - 16 e^p^z^ + 16 p+e^z'' + 60 p+zV- 

+ 24 e^zVi - 88 e^zV+P- + 96 ep| + 104 e^z^p^pS^ + 34 ^2 _ ^2^2^^ ^3^ _^ 276 p+z^p- 

+ 198 p+z^p^ + 156 p+z'^p- + 504 ez^p\pi - 528 ez*p+P- + 912 ezV+P^ - 456 ez^p\p_ 

+ 16 e^zVi + 72 ez^p\ + 16 p^e^z^- + 84 piz^ + 64 p+e^zV- + 96 piz^ - 216 ez'^p^^p+ 

+ 60 p!.z2p+ - 528 ez4p+ pi - 456 ez^p^_p+ + 48 p+ e^z^l - 36 p+z^^ 

- 112 p+e^zV-+ 96 P+ + 276 piz^p+ + 168 ezVt- (B3) 

At a glance, Eq. HB2() appears to depend on the parameter e; however, we can show that it is independent of e, i.e., 
it is conformally invariant. 

2. Contribution of the conical branes 

We can also evaluate the part coming from the cones. Starting with the normal derivatives 



in the bulk we obtain 



.2 1 1 



a,b 



( - 4zV- + 12 p+ p- z^ + 2p^e2^2 ^ 2e2z2p2^ + 24 p+ p_ z^ - 4 p^z^ 



-4 z2p2^ + 12 p+ p_ - 4 p+ e^z^p_ - 4 p2 



- y (2^ (p, + p_ ,2)2 (4 P.^ + 6 P. P- + P^- + ^V- + 4z Vi) . (B5) 
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These values are invariant for the Hmits, z — > 0, oo. The extrinsic curvature in each direction is given by 

= - ^ i^pf 5(1 + z')d. lni2prS,, , = {B6) 

and we find 



2 _ oJ'^K^^ 



(B7) 



Thus, we see that these terms vanish at the locations of the conical branes: at z = and z = oo. 

Given this fact it is easy to show that 04(7 — 1) is also conformally invariant. The explicit form of i?4(r, (5+) for 
conical branes is finally obtained as 



branes 



2887r 27r 1 _ ^ 



o2 £± 



11 + fl_£+)2 

4-^ ^^ , (l-3r) 

(5. 



£±_(l_^)^+(l_£±)r 1^2^(l_|i)2 



where we have used 



2887r r 1 _ ^ L' 5 d _ £±)2 

^ 27r V-^ 271-'' 



5+ - 27r(l - r) 



(1--) 



(B8) 



(B9) 



We have also confirmed that the bulk part of a6(/ = 1) is conformally invariant and we expect a similar behavior 
for the cone part, as in four dimensions. 



APPENDIX C: DERIVATION OF COCYCLE FUNCTION FOR CONICAL BOUNDARIES 

In this Appendix, we shall give a derivation of the cocycle function for conical boundaries in four dimensions. In 
deriving such a term, we shall use the conformal properties of the heat kernel coefficients, see e.g., • 



Conformal properties of the heat kernel coefficients 



In general, for a conformal transformation of the metric, given by 
d-dimensions satisfies the following relation 



6 9ij I 



the heat kernel coefficient in 



Ad) 



[gij,Suj] 



1 



d-2k 



(CI) 



We shall use this formal relation to derive the smeared heat kernel coefficients. At a glance, there seems to be a pole 
at d = 2A; in Eq. I|C1|I . However, because a^^\f = 1) is conformally invariant in four-dimensions (see Appendix B), 
then the leading order pole term vanishes. If a"^\f — 1) is not conformally invariant, then we should add the term, 
'5'-^]U=0j which includes mass and curvature couplings [4^ . This is in order to keep conformal invariance on 

the right-hand-side. In our case a^if\f ~ 1) is already conformally invariant, because of the E term that appears as 
a result of the conformal transformation of the classical scalar action. 
The conical contribution to asjyj — 1) is given by 
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and in our case the conformal transformation is defined by Eq. ljA13|l . 
For a general d-dimensional spacetime, we obtain 



E4e) = 



d-2 



^ g-2(e-l)c:j 



(e- l)Acj- 
r/d- 2 



(e- l)'V*wV.w 



(^)(e-l)Ac.+ (^) (.-1) 



- l)^V*wV,w 



(C3) 



where in the second step we have used the conformal scaling property of the Laplacian, e.g., see 45j. Note that e = 1 
corresponds to the original metric for the codimension two brane. Given the above conformal transformation the 
Ricci scalar, Ricci and Riemann curvature terms scale as [s^, 

1 



Ed + -RW) 



-2(£ 



)2 - 



i)-[5:i?„„-2^i?.,,, + i^«:(-)«:W 

a a.b a,h 

-(e - 1) ((d - 4) ^ Lo,,jn^'''^'n'^''^' + 2/^2^^ + (d - 2) ^ n^°-'^Lo,,n^'''^'^ 

a a 

l)2|^i(rf_4)(rf + 2)^w.,cjyn('')*n(")^' - 2(d - 4)cj.feCj'' 

a 

-2(e-l)^ ^ J^(a)2 _ 2k^^\'^^^ - 2(d - 4)(e - l)w;,n(") V^'"' 



+ (d - 2){d - 4)(e - l)2tj.,w.j-n('^)*n('')^' 



(C4) 



where ; denotes the covariant derivative with respect to the metric e = 1 and curvature tensors which appear on the 
right-hand-side are also defined for the same metric. Here, A2 represents the Laplace operator defined by 



(C5) 



Note that terms including A2W and k'^°'^ lj ■^irS°'^^ vanish on the cone in the model which we are considering and thus, 
we neglect these terms. All the remaining terms proportional to (e — 1) are also proportional to (d — 4). The terms 
which include (or k^")) are also vanishing on the cones. However, in order to see the conformal properties of 

these terms, in this subsection, we keep them. 

By multiplying the conformal transformation of the trace of the brane induced metric, ^/h(e) := e(^-2)(e-i)c^^^ 
with these combinations of intrinsic and extrinsic curvature terms, we can expand the terms as a series in powers of 
((d — 4)). For instance, expanding the first term in Eq. (|C2|I we obtain 

^)\Ea + )r + \,Y^ - 24^)4 
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Similarly, expanding the second piece in Eq. (jC2|) leads to: 

^ - 2 ^ Rata, + i E + ^2 E {^^''^' - 

a a.b a a 
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0((d-4)2). 



(C7) 



2. Cocycle function for the conical branes 

Now we can construct the cocycle function for the conical boundaries. The \/hLLi.jUj'^ terms can be reduced to 
J2a^;i^d^^'^^^^^'^^'' ^^"^ because 



V (p++P-02)2 



vanishes on the cone, at z ^ 0, oo, we can discard these terms. Whence 
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Finally, performing the integration over e we obtain the cocycle contribution for the conical branes: 
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as 



3. Total cocycle contribution 

Combining the above result with the bulk piece, we obtain the total cocycle part of ^4(r, (5+), defined in Eq. (|4.3|1 . 



^4, cocycle (f; <5+) 
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(l + z2)(l+rz2) (l + rz2)4 



21 



9' '^+2-^ 
1447r 27r 1 _ £± 
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1447r 
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where ^'4(1, r, e, z) is given in Eq. I|B3|I . 



(Cll) 



APPENDIX D: WKB ANALYSIS IN THE CONFORMAL FRAME 
1. WKB mass spectrum 

In this Appendix, we derive the mass spectrum in the conformal frame (e = 0): 



and we consider the corresponding eigenvalue problem 



-de^'+dx^ 



(A4 + ^(0)) = -A2(/.a , 
where E{0) is obtained from Eq. H3.16(l and we shall decompose the mass eigenfunction as 

{2n) 



(Dl) 
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(D3) 



From now on, we concentrate on the radial mode function and derive the mass spectrum in the context of the WKB 
approximation. Changing variables to 



2\ 1/2 



Z(l + z2) 

leads to the diagonalized equation of motion in terms of /: 

f"{z) + q{z)f{z)=0 

where 



q{z) = 



(1 + rz2)2Q,2„2 
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2tt 
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(D4) 
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To apply the WKB method we require a Langer transformation ^46J y = \nz and a further redefinition of the 
eigenfunction F — e^^^^f. Thus, Eq. ljD5|) becomes 



{dl+Qiy))Fiy) = 0, 



where 



Qiy) 



A2 - fc2 



g^{ey + e-y)^ {ey + e'y)^ 

The coordinate y runs from —00 to cx), while the turning points, defined by Q{y±) = 0, are given by 



(D7) 
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where y± satisfies ?/_ > > (In the original frame, e = 1, the turning point problem is even more difficult to 
solve.) Furthermore, for the n = modes the turning points extend to y —^ ±oo and in such a case the mass spectrum 
becomes exact in the context of the WKB method. For n ^ 0, turning points appear if and only if 



(DIO) 



The mass spectrum, which will be determined later via the WKB method, is certainly valid for the modes that satisfy 
the above condition. However, we will use the WKB result for deriving the full zeta function, because this is the best 
we can do at present. Nevertheless, it should give a useful order of magnitude approximation to C'(0, ^4:,e=o)- 
In the WKB approximation, the mass spectrum is given by the quantization condition 



1, 



y+ 



(Dll) 



where the integration can be performed as follows: First, we change the integration variable 
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where u = e^^. Then, we use the fact that 
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In the case that we are now considering, 
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{2a -b) + {b- 2c)w_ 
(1 + u-)^b'^ - Aac 



— arcsm 



(2a -b) + {b- 2c)u+ 
(1 + u+)\/b'^ - Aac 



.(D17) 



Unfortunately, in order to find the approximated mass spectrum, we must make one further approximation: A 3> 
max(n, fc). This limit corresponds to the case where (6,m+) — > 00 and u_ — > respectively. Hence, 



\/Q{y)dy 



y+ 



{1 + r)TTan -\ -v/ A^ — /c^ + g^a^ (1 — r^v?- = (to H — )7r . 

g ■ J 2 



(D18) 



The modes with n = correspond to the limits y± — > ^00 and a, c — > 0. Thus, from Eq. I)D17(I the mass spectrum in 
the WKB approximation becomes exact: 



}? = fc2 + g2(2TO+ 1)2 



(D19) 



For these modes we use a single summation in terms of to. 
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2. KK limit 

To investigate the accuracy of WKB method we can compare the result to the non- warped KK hmit, where r — I. 
This has an exact analytic solution for the mode functions with 

= API: if) + BP-^^{f) ""'^TTf^ (°20) 

where we have made a coordinate transformation from z — > f and 

1^(1^ + 1) = ^ (fc2 _ A^) with fi^an. (D21) 

The exact mass spectrum is then determined by the regularities of the mode functions on both boundaries, at f = ±1 
and thus, we find 

I + II + 1/ ~ —771 or 11 — ly — —777 , (D22) 

where 777 = 0, 1, 2, . . . as before. In order to have positivity of the eigenvalue. A, we shall take the latter condition. 
From this the exact mass spectrum is found to be 

A^ = g'^{2m + 2an){2m + 2 + 2an) + k'^ 

= g^{2m + 2an+lf + - . (D23) 

The result of our WKB method, Eq. (|D18|I . agrees with the one above (when r = 1) except for a constant term in the 
eigenvalue spectrum. However, given that our WKB approximation is valid for large quantum numbers this constant 
term is negligible in such a limit. 

3. WKB zeta function 

Given the WKB spectrum derived in the previous subsection, we derive an expression for the zeta function and the 
associated effective action. The zeta function is given by the mode sum 

(27r)2CwKB(s,A4,.=o) = j d^^Yl f d^kX'^" ■ (D24) 

Hereafter, we omit "WKB" in the subscript of the zeta function. Here, we shall use the density of states method [4^ . 
which from Eq. (|D18p the density of states is easily derived: 

- ^ . ^ • (D25) 

dX 2g ^A2-/c2+772g2Q,2(l_r)2 

We can therefore replace the sum over m with an integration in terms of A ^3 

m 

It the following it will be convenient to decompose the zeta function into two pieces, i.e., the 77 7^ and 77 = modes. 
The total of the derivative of zeta function is finally given by 

C'(0, A4,e=0) - Cn#o'(0, A4,e=o) + Cn=o'(0, A4..=o) • (D27) 

Similar considerations can easily be applied to the six-dimensional case, working in the conformal frame. 
For the 77 7^ mode, we obtain the following expression 



d^x V ' / dXX / dkk{2^)( — )x 



n— — 00 
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where because of the condition that there be two turning points, see Eq. (jD10|) . we must restrict the range of 
integration. The ' denotes that the n — mode is to be dropped, because the potential for this mode is not oscillator- 
Uke. Performing the A:-integration first, we obtain 

n=l ^ 

OO nOC 

- 2^(-) / d~X~x(^X'' +4ra^g^n^y\ag{l + r)n)} , (D29) 

n=l ^ 

where A is defined by 

= - Ara^g^n^ . (D30) 
For sufficiently large Re(s), we obtain the following integration formulae 



dAA A^ +4m2g^n^ 



-J—^iArr-^{agnf'-'\ (D31) 



dXx(x^ + 4ra^g^n^^ " ^X^ + a^g^{l + r)^n^ 

where again 2Fi{a,b,c; z) is a hypergeometric function. Note that the first term in Eq. (jD32|l has been analytically 
continued by employing the standard duplication formula for the Gamma function; otherwise the result will be ill 
defined when taking the derivative at s = 0. 

Finally, expressing the summations over n in terms Riemann zeta functions we obtain (substituting back the 
expression for a, see Eq. (C6)) 



\3-2s 



(4r)i-«(l + r)/ ^ , 1 , , ^. 4r 



2(s- 1) V" 2' ' ' (l+r)2 

where C-rX'S) is Riemann's zeta function. Then, taking the derivative in terms of s and an analytic continuation s — s- 0, 
we find 



C#o(0,A4,.o) = 270(2/- ^+)3 / d^{-2-3j-m-l)){l-rr + 2(-3 + ry{-3 + 3ln{^^)) 

4r 

(TT7) 



+ (1 + r)(l - 4r + r^) (e In (^g^^) " 720C;,(-3)) + 18(1 + r)r-^^2F^ [-1,1, 2, 7Y^^34) 



The zeta function for n = mode is given by 

{27TfCn=,{s, A4^,.o) = fd'xf <fk [5^(2™ + 1)2 + = ^(M!^(^^(2s - 2, i) , (D35) 

m=0 ^'^ 

where ^//(s,a) is Hurwitz's zeta function. Taking the derivative with respect to s and continuing to s — > gives 

W(0,A4.e=o)-^^^%^. (D36) 
The total resuh of Eq. (g^l or Eq. (jO^ is 

^4,WKB(r, <5+) = - r^^Q^fjls^^s {{ - 2 - 37 ^ 3V'(-i)) (1 - rf + 2(-3 + ( - 3 + 3 In ( (Y^)) (^37) 

9^/ /27rg(l-r), , . ^^ . ,9,1 4r ,^ 3g'^Cn'{-'2) 

+ (1 + .)(1 - 4. + r^) (6 In {^^) - 720Ck(-3)) + 18(1 + r)r-,F,[--, 1, 2, ^^]} + . 
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